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Abstract 

Using data collected by the FOCUS experiment at Fermilab, we report the discovery 
of the decay modes D° — > K~iT + iT + Tr + Tr~iT~ and D° — > 7r + 7r + 7r + 7r~7r~7r~. With a 
sample of 48 ± 10 reconstructed D° — > i"C - 7r + 7r + 7r + 7r - 7r - decays and 149 ± 17 
reconstructed D° — > 7r + 7r + 7r + 7r _ 7r~7r~ decays, we measure the following relative 
branching ratios: 

T{D° -» ^-vr+7r+7r+7r-vr-)/r(L> -► K-7r+7r+^-) = (2.70 ± 0.58 ± 0.38) x 10~ 3 

T(D° -» 7r+7r+7r+7r-7r-7r-)/r(L> -» ^-vr+Tr+vr-) = (5.23 ± 0.59 ± 1.35) x lO" 3 

r(D° -» ir+iT + iT+iT-iT-iT-)/r{D -» ^-vr+Tr+Tr+vr-vr-) = 1.93 ± 0.47 ± 0.48 

The first errors are statistical and the second are systematic. The branching fraction 
of the Cabibbo suppressed six-body decay mode is measured to be a factor of two 
higher than the branching fraction of the Cabibbo favored six-body decay mode. 



1 Introduction 

Hadronic decays of charmed mesons have been extensively studied in recent 
years. However, six-body hadronic decays of the D° have not been previously 
observed; only an upper limit exists for the D° — > 7r + 7r + 7r + 7r~7r~7r~ branching 
fraction [1]. In this paper, we present the first branching ratio measurements 
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of the D° — > A' _ 7r + 7r + 7r + 7r _ 7r _ and D° — > 7r + 7r + 7r + 7r~7r~7r~ decay modes. 
Charge conjugate states are implicitly included and we use the abbreviations 
D° — > K5n, D° — > 67r, and _D° — > i^37r for the fully charged states. 

The fixed-target charm photoproduction experiment FOCUS collected data 
during the 1996-1997 fixed-target run at Fermilab. The FOCUS detector is a 
large aperture spectrometer with excellent vertexing and particle identification 
capabilities. A photon beam is derived from the bremsstrahlung of secondary 
electrons and positrons produced from the 800 GeV/c Tevatron proton beam. 
The photon beam interacts with a segmented beryllium-oxide target. The av- 
erage photon energy for the interactions collected for the measurements we 
report is 180 GeV. Charged particles are tracked by two systems of silicon 
microvertex detectors. The upstream system [2], consisting of 4 planes (two 
views in two stations), is interleaved with the experimental target, while the 
other system lies downstream of the target and consists of twelve planes of 
microstrips arranged in four stations of three views. These detectors provide 
high resolution separation of production and decay vertices. The momentum 
of a charged particle is determined by measuring its deflections in two anal- 
ysis magnets of opposite polarity with five stations of multiwire proportional 
chambers. Three multicell threshold Cerenkov counters are used to discrimi- 
nate between electrons, pions, kaons, and protons. 



2 Signals and selection criteria 

A candidate driven vertexing algorithm [3] is used to reconstruct D° decays 
into six-body final states. A D° candidate consists of six tracks in an event 
that have zero total charge and form a vertex with at least a 2% confidence 
level. The momentum vector of the D° candidate, formed from the momenta 
of the six tracks, is then intersected with at least one other track in the event 
to form the D° production vertex; the confidence level of this vertex is re- 
quired to be at least 1%. Additional cuts are applied based on event geometry 
and particle identification. To minimize systematic errors, identical cuts are 
used on the two six-body decay modes and on the normalizing mode, except 
that the D° — > 6tt mode has no kaon identification cut. Our most effective 
cut for reducing non-charm backgrounds is a significance of detachment cut 
that requires the separation, £, between the D° production and decay ver- 
tices divided by its error, a^, to be greater than some threshold, in our case 
Ijai > 13. The D° decay vertex is also required to be located outside of ma- 
terial in the target region by at least 2 standard deviations, which serves to 
reduce backgrounds from secondary interactions. Cerenkov particle identifica- 
tion is done using a x 2 -like variable Wi = — 2 lnLikelihood(i), where % ranges 
over electron, pion, kaon, and proton hypotheses [4]. For each pion candidate, 
we require min{lU e , W K , W p } — W w > —4, which requires that each pion can- 
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Fig. 1. Invariant mass distributions of (a) K5tt and (b) 6tt. 

didate is not highly favored to be an electron, kaon, or proton rather than a 
pion. For the kaon candidate in D° — > K5n and in the normalizing mode, we 
require that the kaon hypothesis is more likely than the pion hypothesis with 
the cut W n — Wk > 3. Finally, the largest confidence level that one of the 
tracks from the decay vertex intersects the production vertex is required to be 
less than 25%. 



The invariant mass distributions of the D° candidates that satisfy these cri- 
teria are plotted in Figure 1. The D° — > K5n mass plot is fit with a linear 
polynomial plus two Gaussians with the same mean but different widths. We 
fit with two Gaussians because the mass resolution varies with momentum and 
the position of the decay vertex; the sum of two Gaussians provides a much 
better approximation to this situation than a single Gaussian. The widths 
of the two Gaussians and their relative yields are fixed to values obtained 
from a Monte Carlo simulation (61% of the total yield is in a Gaussian shape 
with a = 5.9 MeV/c 2 and 39% of the total yield is in a Gaussian shape with 
a = 13.1 MeV/c 2 ). The fit returns a signal yield of 48 ± 10 events. Based on 
studies of reflections above and below the signal, we choose to fit over the 
range 1.78 GeV/c 2 to 1.98 GeV/c 2 . The reflection below 1.78 GeV/c 2 is con- 
sistent with partial reconstruction of seven-body final states from D + — > KQir 
and from the decay chain D° — > K~ir + r]', rj — > ir + iT~ri, rj — > it + it~it or 
r) — > 7r + 7r~7, which yields the same final state as K57i with an additional 7r° 
or 7. The structure above 2.0 GeV/c 2 is due to D + — > K~7t + 7t + 7t + tt~ with a 
random 7r~ intersecting the decay vertex. 

The D° — > 6n mass plot is also fit with a linear polynomial plus two Gaussians 
with the same mean whose widths and relative yield are fixed to values from 
a Monte Carlo simulation (60% of the total yield is in a Gaussian shape with 
a = 7.6 MeV/c 2 and 40% of the total yield is in a Gaussian shape with 
a = 16.2 MeV/c 2 ). The fit returns 149 ± 17 events. 
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Table 1 

Branching ratio measurements. The first error is statistical and the second is sys- 
tematic. 



Decay Mode 




Branching Ratio 


r(D -*A"-7r+7r+7r+7r- 


-7T-) 


(2.70 ± 0.58 ± 0.38) x 10~ 3 
(5.23 ± 0.59 ± 1.35) x 10~ 3 
1.93 ± 0.47 ± 0.48 


r(D u ^.ft:-7r+7r+7r 
r(D -»7r+7r+7r+7r-7r- 


-) 
TP") 


r(D'J-*if-7r+7r+7r- 
r(D°^7r+7T+7r+7r-7r- 


") 

-7T-) 


r(D°^X-7r+7r+7r+7r- 


-7T-) 



We measure the branching ratios of D° — > K~it + it + it + it~it~ and D° — > 
7T + 7r + 7r + 7r _ 7r _ 7r~ relative to the high statistics mode D° — > 7^ _ 7r + 7r + 7r _ . The 
D° — > fT37r normalizing mode is fit in the same way as the two six-body 
modes, and the fit returns 70 466 ± 277 signal events. We also directly mea- 
sure the relative branching ratio of the two six-body decay modes T(D° — > 
7r + 7T + 7r + 7r~7r~7r~)/T(.D — > K~ty + 'k + 'k + 'k~'k~) in order to take into account 
any correlations in systematic errors on the two modes. From Monte Carlo sim- 
ulations using a nonresonant model of each six-body decay mode, we compute 
the relative efficiencies fjjfefffg = 0.254 ± 0.004, jfefgj = 0.405 ± 0.004, 

an< ^ e(D a ^K^) = 1-^96 ± 0.027. The resulting branching ratio measurements 
are shown in Table 1. Since the D° — > 6n mode is Cabibbo suppressed while 
the D° — > K5tt mode is Cabibbo favored, one might expect the relative 
branching ratio of these two modes to be about tan 2 6 C ~ 0.05, where 6 C 
is the Cabibbo angle. Our measurement of their relative branching ratio is 

t^SkII) = L93 ± °- 47 ( stat -) ± °- 48 ( § y s -)- 



3 Systematic errors 



To check for bias in our selection of cuts, we have studied the sensitivity of 
the results to cut selection by individually varying each cut over a reasonable 
range of values. The branching ratio measurements are stable as the cuts are 
varied. We have investigated a number of sources of systematic uncertainty 
in the branching ratio measurements. These sources are described below, and 
the systematic errors associated with them are listed in Table 2. 

We quantify the systematic uncertainty on the efficiency due to Monte Carlo 
simulation using the split sample procedure based on the S'-factor method 
employed by the Particle Data Group [5,6]. The data set is split into indepen- 
dent subsamples by reconstructed D° momentum and by early and late runs, 
which have different target and silicon microvertex detector configurations. 
Because of our limited statistics, the splits are done one variable at a time. 
We measure the branching ratio for each independent subsample and assess 
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Table 2 

Systematic error contributions as a percentage of the branching ratio. 



Source 


r(D°- 
r(D«- 


>K5n) 
>K3n) 


1 (V -^bw) 
r(D°- > X37r) 


TV nO 


->b7Tj 

>K5tv) 


Run period split 


13.1 


Wo 


25.3% 


24.: 


2% 


Momentum split 


0.0 


% 


0.0% 


o.c 


1% 


Fitting 


3.4% 


4.1% 


5.3 


;% 


Subresonances 


2.8 


% 


2.5% 


3.6 




6-body vs. 4-body 


2.8 


% 


2.8% 






Total 


14.1 


Wo 


25.9% 


25.1 


3% 



whether the subsample measurements are consistent with a single value by 
examining the \ 2 ■ If X 2 / (degrees of freedom) > 1, we scale up the errors such 
that x 2 /d.o.f. = 1. If the scaled error on the weighted average of the subsample 
measurements exceeds the statistical error on the whole sample measurement, 
we define the split sample systematic error to be the difference in quadrature 
between the scaled error and the statistical error. 

We have studied the dependence of the results on the fitting procedure by 
fitting the histograms in several different ways: using one Gaussian instead 
of two, different bin sizes, and different background parameterizations. In the 
D° — > K5n case, we also used fit functions that included reflection shapes 
obtained from Monte Carlo simulations. The sample standard deviation of 
the branching ratio measurements from the different fit variants is taken as 
the fit variant contribution to the systematic error. 

Because the resonance substructures of D° — > K5n and D° — > 6n are un- 
known and the subresonance model used in the Monte Carlo affects the re- 
construction efficiency, we compute the efficiencies for several subresonance 
models and use the sample standard deviation of the resulting branching ra- 
tios as a contribution to the systematic error. The subresonance models used 
for D° -> K5ir are D° -> i^a 1 (1260) + , D° -> K*° p(1450) (p(1450) -> 4tt), 
and D° — > K 4ir in addition to a nonresonant model. For D° — > Git, the 
subresonance model D° — > 7r~ai(1260) + is used in addition to a nonreso- 
nant model. For subresonance models of both modes involving the ai(1260) + , 
we compute the efficiency for each of three different models for the decay 
of the ai(1260) + : ai(1260)+ -> / 2 (1270)7r+, ai(1260)+ -> /o(1370)tt + , and 
a!(1260)+ -> p(1450)tt+, where the / 2 (1270), / (1370), and p(1450) decay to 
four charged pions. The mass and width of the ai(1260) + are assumed to be 
1230 MeV/c 2 and 400 MeV/c 2 , respectively. 

We also include a systematic error contribution from differences in absolute 
tracking efficiencies for six-body versus four-body final states. The total sys- 
tematic error is obtained by adding the different contributions in quadrature. 
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4 Conclusion 



We have presented the first measurements of six-body hadronic decays of the 
D°. The results are: 

(2.70 ± 0.58 ± 0.38) x 10~ 3 
(5.23 ± 0.59 ± 1.35) x 10~ 3 
1.93 ±0.47 ±0.48 



t(d°^k-tt' 


TT+TT + TT' 


TT' 


") _ 


t(d°~>k- 


'TT+TT + TT' 


~) 




r(D°^7r+7r+ 


7T+TT~7T~ 


' 7T~ 


) _ 


r(D°^_fsr- 


TT+TT + TT' 


") 




r(D°^7r+7r+ 


TT + TT~ TT~ 


" 7T~ 


") _ 


T(D°^K-n J 


'TT+TT + TT' 


~7T" 





The relative branching ratio of the two six-body decay modes is much higher 
than one might expect from Cabibbo suppression. Theoretical discussion of 
many-body charm decays has suggested a "vector-dominance model" in which 
a charmed meson emits a W ± which hadronizes into a charged vector, axial- 
vector, or pseudoscalar meson [7]. Studies of five-body charm decays by FO- 
CUS have provided evidence for this model with five-body decays of the 
D°, D + , and Df being dominated by quasi-two-body decays involving the 
ai(1260) ± [8,9]. Our result for T(D° -> Qtv)/T(D° -> K5tt) may be qualita- 
tively explained by the hypothesis that six-body decays of the D° proceed 
primarily through quasi-two-body decays involving an ai(1260) + . The decay 
channels of the ai(1260) + that can result in five charged pions are /2(1270)7r + , 
/o(1370)tt + , and p(1450)tt+. If the D° decays to K~ ai (l260) + , then only frac- 
tions of the widths of the ^(1270), /o(1370), and p(1450) are available for the 
decay of the ai(1260) + , resulting in a significant suppression of six-body final 
states involving a kaon compared with six pion final states from the decay 
D° -> 7r-a 1 (1260) + . 
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National Science Foundation, the U. S. Department of Energy, the Italian Is- 
tituto Nazionale di Fisica Nucleare and Ministero della Istruzione Universita 
e Ricerca, the Brazilian Conselho Nacional de Desenvolvimento Cientifico e 
Tecnologico, CONACyT-Mexico, and the Korea Research Foundation of the 
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